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Abstract 

We present an alternative organizational scheme for developing effective theories of 2- and 3- 
body systems that is systematic, accurate, and efficient with controlled errors. To illustrate our 
approach we consider the bound state and scattering properties of the 2- and 3-nucleon systems. 
0\ Our approach combines the computational benefits of using separable potentials with the improved 

^ convergence properties of potentials evolved with a renormalization group procedure. Long ago 

Harms showed that any potential can be expanded in a series of separable terms, but this fact is 
only useful if the expansion can be truncated at low order. The separable expansion provides an 
attractive organizational scheme that incorporates the two body bound state in the leading term 
while allowing for systematic corrections thereafter. We show that when applied to a renormaliza- 
tion group-evolved potential, the separable expansion converges rapidly, with accurate results for 
both 2- and 3-body scattering processes using only two separable terms. 

PACS numbers: 21.45.+v,34.50.-s,03.75.Nt 
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I. INTRODUCTION 



Effective Theories (ET's) provide a potent and practical approach to understanding low- 
energy nuclear few-body systems where it is assumed that the low-energy long distance 
features should not depend on details of the short distance physics [H El El 111 El E] • Most 
contemporary implementations are based on a field theory specified by a Lagrangian con- 
taining a series of contact interactions with increasing powers of momentum consistent with 
the symmetries of the underlying theory. This series, it is hoped, may be truncated at some 
low order and still account for the low energy physics once the (few) parameters of the ET 
are fixed by external inputs. The power of the approach is that these few inputs yield a 
theory which can then be used to calculate systematically any low energy property of the 
system with relative ease. 

However, the use of contact interactions as a starting point, while nearly ubiquitous, is 
not essential. The original motivation for this choice was presumably one of simplicity. If the 
physics to be described were perturbative in nature, the organization of the theory would be 
trivial as terms in the series of contact interactions would appear at the order of their naive 
scaling dimension. However, the presence of a weakly bound 2-body state in the ^5*1 channel 
- the deuteron - renders the problem fundamentally non-perturbative and the organization 
of contributions from contact terms in the ET's is, as some practitioners acknowledge, far 
from simple and often controversial. Also, starting with contact interactions appears to be 
awkward if one wishes to treat finite range effects in the 2-body interaction. Indeed, Phillips, 
Beane, and Cohen [7] have shown that a 2-body potential consisting of contact terms and 
derivatives thereof necessarily leads to a negative range parameter in the Effective Range 
Expansion (ERE) of pcot{S) regardless of the number of terms in the expansion. Such 
considerations lead one to consider alternative organizational principles that acknowledge 
the finite range of the 2-body interactions at the outset. One such approach is based on 
treating the 2-body interaction as a sum of separable terms, each consisting of the outer 
product of form factors which reflect the structure of real or virtual 2-body states. 

Separable potentials have a long history in nuclear physics going back to the period 
immediately following Faddeev's seminal work on the 3-body problem when calculations 
of Nd scattering relied on finite-range separable A^A^ interactions to an extent that was 
as nearly universal as the use of contact interactions in ET's is today. More recently, 
Lepage [8] provided a general approach to developing effective interactions starting with 
the now-canonical series of contact interactions but morphing into a similar series smeared 
in configuration space by what amounts to a form factor. The possibility of using separable 
potentials is also mentioned in passing in several of the important papers in the recent 
history of ET's addressing three-body physics including Afnan and Phillips |9] . In any 
case, there is much current interest in how best to treat finite range effects that may enter 
in Efimov physics which is now accessible experimentally in cold atom systems where one 
exploits Feshbach resonances to tune the 2-body scattering length, a, to infinity. In this 
regime it would seem prudent to use an approach that does not implicitly assume range 
effects are negligible in leading order. 

To establish the methodology in a simple context we use the uncoupled ^Si potential 
of Malfiiet and Tjon (MT-III) [10] which enables us to treat the 3-nucleon quartet (^5*3) 
channel, while avoiding complications due to tensor coupled phase shifts and Efimov physics. 
Extension to the doublet case will be treated in a separate paper. In this work, we exploit 
the work of Harms [TT] who showed that a general 2-body potential can be written as an 
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infinite sum of separable terms. Of course, the utility of the approach depends critically 
on the convergence properties of the expansion, a study of which is the focus of this work. 
There exist fully general Faddeev treatments of the nd system based on realistic microscopic 
potentials, see [12], [T3] and references therein. The goal of this work is to provide a com- 
putationally simpler yet accurate alternative approach exploiting both separable expansions 
and renormalization group methods. 

We consider here three approaches for implementing the Harms expansion. First, we 
directly expand the bare 2-body potential to obtain the separable representation of the 
potential. We find the convergence of the separable expansion for the bare potential is 
reasonable, but only if a judicious, but somewhat arbitrary, ordering of the separable terms 
is made. In addition the Harms form factors for the bare interaction have considerable 
support at large momenta requiring a large cut-off in the loop integral. Problems with non- 
physical negative imaginary phase shifts arise when the loop integration cut-off is reduced. 
This suggests that one may profit by first integrating out the high momentum degrees 
of freedom using a renormalization group scheme, such as the Similarity Renormalization 
Group (SRG). 

The use of the RG in nuclear calculations has been extensively explored in a series of 
papers by Bogner and collaborators [HI [151 US [H!- More specifically, Bogner, Furnstahl, 
and Perry [16] have recently shown how the SRG drives the matrix elements of the bare 
A^A^ interaction to the diagonal and, in a subsequent paper, used the SRG in the analysis 
of a 2-state model 3-body problem [T7]. Thus, our second approach is to first evolve the 
bare A^A^-potential to some low momentum scale before expanding the resulting Viow using 
the methods of Harms. We find that the convergence properties of the separable expansion 
are much improved, requiring only two terms to reproduce the exact 2-body T-matrix out 
to about 400 MeV/c and suffering no issues with negative imaginary phase shifts in 3-body 
applications. 

Our third approach is to develop a phenomenological separable potential (dipole or gaus- 
sian times polynomial form) and fix its parameters by fitting to the 2-body phase shifts. In 
an earlier work [? ], we used this approach by truncating the separable expansion at lead- 
ing order (an approximation historically referred to as the "Unitary Pole Approximation" 
). We then used a phenomenological dipole form factor with parameters fit to the phase 
shifts. By expanding this model in the inverse of the dipole cut-off parameter, we were able 
to study systematically finite range expansions, to make connection to the modern ERE's 
based on contact interactions, and thereby to reveal limitations of the ERE in accounting 
for finite range effects in the Nd system. Kamada, et al. also studied 3-nucleon systems 
using separable expansions of low-momentum effective potentials using a modified Suzuki 
and Okubo approach to generate Viowi but using a phenomenological separable expansion in 
Legendre polynomials [I^. Their expansion converged relatively slowly; requiring as many 
as 16 terms in the separable expansion. 

We next move to an examination of the 3-nucleon problem, confining our attention to 
the relatively simple quartet, ^5*3, case to develop our methodology. We confirm the utility 
of combining the SRG and the Harm's separable expansion by calculating three nucleon 
quartet phase shifts, including the sensitive and small imaginary phase shifts, accurately 
using only rank-2 approximation. We emphasize that the approximation method is general, 
systematic, and controlled. 

The paper is organized as follows. In Section II we review Harms' derivation of the 
separable expansion which defines our notation. We apply the method to the Malfiiet-Tjon 
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(MT-III) potential, Ref- ^Q\, and show its convergence properties as the rank of the expansion 
is increased. In Section III we review the Similarity Renormalization Group and then use 
it to flow the bare MT-III NN potential to low momentum. In doing so, as also shown 
by Bogner et al., we illustrate how the SRG suppresses the off-diagonal components of the 
potential at high momentum thereby effectively rendering the high momentum components 
of Viow sterile to scattering. We conclude Section III by expanding the SRG-evolved potential 
in a separable expansion and show that the convergence properties are improved over that 
of the bare potential. In Section IV we turn to the 3-body problem by treating Nd quartet 
scattering. We first review our UPA results and then examine those of the rank-2 separable 
approximation to the SRG evolved MT-III potential. The agreement with the full Faddeev 
calculations of Huber, et al. [12] and Kievsky, et al. |T3] is excellent. We discuss our results 
and prospective future developments in the concluding section. 

II. HARMS' SEPARABLE EXPANSION 

In the early days of nuclear 3-body calculations, practitioners made use of a simple 
separable expansion of the 2-body interaction to render the 2-body problem in anaytic 
form [ini EHl EH [221 ESI EH ES]. These applications employed the simplest form of a 
separable interaction, namely one consisting of a single term which, in terminology later, is 
a rank-1 separable potential. These early papers emphasized that such an approximation is 
well-suited to situations where the 2-body properties are dominated by a single real or virtual 
bound state, or dimer. Use of a rank-1 interaction was typically referred to as the Unitary 
Pole Approximation (UPA) since, as will be shown below, 2-body unitarity is guaranteed 
using this form of the interaction. 

In his seminal 1963 paper, Weinberg [19] showed how the eigenvalues of the operator 
Go{E)V can be used as an organizing principle for developing a perturbation expansion of 
the T-matrix. In 1970, using this organizing principle, E. Harms showed how any 2-body 
interaction can be expressed exactly as an infinite sum of separable terms. Truncation of the 
expansion at some finite rank allows the construction of the fully off-shell 2-body T-matrix by 
inverting a matrix with a dimension of the rank retained. His expansion yields the UPA when 
a single term (associated with the bound state) is retained. It should be noted that since 
separable expansions are always compact, they cannot converge for local potentials which 
are non-compact [261EI]- Nevertheless, as pointed out by Koike [2H], a separable expansion 
may approximate local potential quite well in limited regions of momentum and energy, and 
specifically for applications to 3-body problems the Green's function in the connected-kernel 
3-body equation regulates the integrand rendering the separable expansion convergent for 
this important case. We now briefly reconstruct Harms' arguments using his notation with 
one exception specified below. 

Consider the 2-body Schroedinger equation (^, c = 1, nineutron = f^proton = f^)'- 



where Gq{E) is the free 2-body propagator. Assume there exists a single bound state at 
E —B = — The bound state satisfies 



G^\E^)\Xn) = V\Xn) 



(1) 




(2) 



|Vn)=Go^(-i?)|Xn). 



(3) 
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We see that the bound state Schroedinger equation can be written as 

VG^{-B)\i,B) = \B\iJB) (4) 
where Xb = 1- Solutions \il)n) of the generahzation of Eq. |4] 

VGQ{-B)\i;n) = \n\i^n) (5) 

form a complete set of kets with normalization condition 

mG^{-B)\i}n) = -Km- (6) 




FIG. 1: (Color online). pcoi{5) in the channel for potential MT-III. The black diamonds are 
the full Lippmann-Schwinger (LS) calculation taken to be exact. The light dotted (red), short 
dashed (red), and light (red) lines are the rank-1 (UPA), rank-2, rank-4 results, respectively, for 
the Harms separable expansion of the bare potential, Vq- The heavy dotted (blue), long dashed 
(blue), and solid (blue) lines lines are the rank-1 (UPA), rank-2, and rank-4 results, respectively, 
for the Harms separable expansion of the SRG-evolved potential, Vio^. 

We call the |V'n)'s the "form factors" of the potential V . Note that our A„'s are the inverses 
of those defined by Harms. Also note that, with this normalization convention, the bound 
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state wave function is not normalized. We may define the wavefunction renormalization Zb 

{4jb\Go{E)\^b). (7) 



via 

Zb' = (XbIxb) = iiJBlGoi-BYl^B) = - 



Wfien there is a single bound state, the eigenvalues satisfy 

An^B < Xb = I- (8) 

Weinberg noted that the eigenvalues can be used to develop a perturbative expansion of the 
T-matrix which will converge when |A| < 1 [19]. However, for the MT-III potential there 
are large negative eigenvalues; so the convergence of an expansion based on this ordering 
principle is in question. We observe that the potential can be written as 

n 

which is verified by substitution into Eq. |5] and then using the normalization condition |6} 
The 2-body T-matrix satisfies the Lippmann-Schwinger (LS) equation 

T{p,p;E-) = V{p,p) + J ^ E^-qym ^^^^ 
It is easy to show that 

T(p, p'; E^) = -J2 (Pl^n) A(i?+)„,^ (^^IpO, (11) 

n,m 

where 

[A-\E^)]n,m = \-'Sn,m + |Go (^+) |^„.) • (12) 

Thus, we see that the 2-body T-matrix associated with a separable potential to a given rank 
involves inverting a matrix with dimension of that rank. The UPA corresponds to the single 
term, 

VuPA ^ -I^b) O^bI (13) 

where Ab = 1 has been used, and 



Tupa{p,p';E^) -(pI^b) Aupa{E^) {^PbIp') = "^^'^Vr^lJ^ . (14) 



(pI^b) ii^B 


p') 


i + ii^B 


Go{E+)\ijB) 



which we recognize as a standard result for (rank-1) separable 2-body interactions |29] . 
The residue of the T-matrix at the deuteron pole, p — 27, is related to the wavefunction 
normalization, Eq. [7j 

'JZ[TuPA]\p^r,^\{il\i^B)\'ZB. (15) 

As with our predecessors, our interest in separable forms of 2-body interactions stems 
from the great simplification they afford in solving 3-body Faddeev equations ([29]). In 
the context of the Harms expansion, when focussing on low energy properties of 3-body 
systems, the promised simplicity is only realized if the separable expressions for V and T 
given above can be truncated at some low order and still be accurate over the relevant range 
of momentum and energy. We now examine this issue using the familiar ^Si NN potential 
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Potential 


A > 


A < 




1.0000 
0.159190 
0.061669 
0.032467 


-1.8917 
-0.24285 
-0.018806 


Vlow 


1.0000 
0.11937 
0.024633 
0.0055180 


-0.11781 
-0.05982 
-0.050676 



TABLE I: First few positive and negative Harms eigenvalues for Vq and V^oui for the MT-III 
potential. 



of Malfiiet and Tjon [TO]. We use this potential because we consider the present work to 
be a "proof of principle" exercise, and we wish to avoid the complexities of the full coupled 
channels nature of the ^5*1 — ^ Di system while treating the quartet channel in the nd system 
where only the NN triplet channel contributes and complexities involving Efimov physics 
are absent. 

In Figure [l| we compare pcot5 for Potential III of Malfiiet and Tjon (MT-III) computed 
by exact solution of the LS equation (Eq. 10) with results of the T-matrix (Eq. 11) using 
the Harms method for various truncations. Ignore for now the curves involving "VJou," . We 
always take the leading term to be that associated with the bound state, Ai — > A^ = 1. We 
order the separable expansion by the eigenvalues, A„, of the operator GqV , Eq. |5| separately 
tracking positive and negative terms, see Table [Tj Weinberg showed this forms the basis for 
a systematic perturbation expansion in that it orders the separable terms according to the 
magnitude of their contribution to the potential [HI [TH]. Figure 1 shows the convergence 
of the separable expansion of the bare potential as the rank is increased, first from rank-1 
(UPA), rank-2, fully converging by rank-4. (By a judicious, but arbitrary, choice of ordering, 
convergence can actually be achieved at rank-3.) We see that for either potential the rank- 
4 Harms results are essentially in perfect agreement with the full LS solution over the 
momentum range shown (p^ < 4 fm~^). (This finding is of some interest in its own right 
as the numerical cost of the Harms calculation is significantly less than that of the full LS 
computation.) 

Figure 2 shows the form factors for the first few terms in the Harms expansion of the bare 
potential. As can be seen, aside from the first (UPA) term, all the higher rank form factors 
have considerable support at large momenta. (Ignore for now the curves involving "VJo^".) 
Truncating at rank-1 means that the calculated phase shift will go to zero too rapidly with 

giving a pcot5 consistently above the exact value as shown in Figure 1. To capture the 
high momentum dependence therefore requires that higher order terms be retained in the 
Harms expansion. As we are interested only in the low momentum physical properties, this 
suggests that a renormalization group procedure could improve the convergence properties 
of the Harms expansion. This possibility is explored in the next section. 
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FIG. 2: (Color online). Form factors, {p\'>pn): from the solution of Eq. |5] for the "^Si channel. 
The light solid (red) curve is the bound state form factor (A = 1.00000) and the short dash (red) 
and short dash-dot (red) curves are the second (A = 0.159190 . . . ) and third (A = —1.891665 . . . ) 
form factors for the Harms separable expansion of the bare potential, Vq. The thick solid (blue) 
curve is the bound state form factor (A = 1.00000), the long dash (blue) curve is the second (A = 
0.1193678 . . . ) form factor, and the long dash-dot (blue) curve is the third (A = 0.024633 . . . ) form 
factor for the Harms separable expansion of the SRG potential, Viow- Note the strong supression 
of the SRG- (blue) evolved form factors for p > A^ou, = 2 fm~^. For comparison we also show the 
dotted (black) curve is the phenomenological dipole form factor. 

III. DETERMINATION OF Viow VIA THE SIMILARITY RENORMALIZATION 
GROUP 

Much recent work has focussed on using a renormalization group approach to develop 
effective 2-body interactions - designated as Viow ^ evolved from "bare" interactions defined 
at an energy scale much higher than the scale of the physical application in which the 2-body 
interaction is to be used. Several methods for obtaining Viow from the bare potentials have 
been developed over the years including the Kuo-Kurasawa (KK), Lee-Suzuki (LSuz), Con- 
tractor Renormalization (CORE) and Similarity Renormalization Group (SRG) approaches. 
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We have performed calculations using each of the last three methods and find similar results 
but here focus on the SRG whose basic features we now review. 

The basic concepts of the SRG are straight-forward; they were formulated independently 
in the early 1990's by Wegner [30] and by Glazek and Wilson [211 132] • One begins with a 
Hamiltonian for the 2-body interaction, Hq = T + Vq. For simplicity, we take a basis of 
momentum eigenstates with a maximum momentum, Aq, assumed to be very large com- 
pared to the momentum scale of the 2-body physics of interest. Finding the eigenvalues 
and eigenvectors of Hq constitutes a full solution of the 2-body problem. The SRG allows 
determination of Hiow = T + Viow via a sequence of similarity transformations which (i) 
retains all original eigenvalues while (ii) reducing the upper limit of momentum for which 
Viow has non-trivial dynamical content - which is identified with Aq for Vq - to some lower 
momentum scale, Aiow- Consider the similarity transformation, 

H{s) = U{s) Ho U{s)\ (16) 

where s is a scale parameter related to momenta A between Aq and Aiow Differentiation of 



Eq. [16] with respect s yields 

^-^ = [vis), His)] (17) 

where 

vis) = ^ U^is) = -vis)'. (18) 



We recognize Eq. 17 as a Renormalization Group (RG) flow equation for His). The matrix 
rjis) specifies the flow and a good choice for r]is) in our momentum basis where the kinetic 
energy operator T is diagonal has been shown to be 

vis) = [T,His)] (19) 

which yields the familiar double-commutator form of the SRG flow equation 

^-^=[[T,His)],His)]. (20) 
Writing out this equation explicitly in terms of matrix elements of T and His), we find 

= —iTii — Tjj) Hij + 2_^l^ii + ^ii ~ '^Tkk\HikHkj (21) 

k 

where we have dropped the explicit dependence of H on s. Focussing on the first term on 
the rhs of this equation, we see that the SRG "flow" drives off-diagonal matrix elements of 
H toward zero at a rate which is quadratic in the distance from the diagonal. This behavior 
dominates even when the contributions of the second term are included. Thus the SRG flow 
(i) keeps all original eigenvalues (i.e., at Aq) unchanged since successive if's are related by 
a similarity transformation and (ii) drives the Hip,p') matrix toward a diagonal form for p 
or p' > Alow. Hence, these high-momentum modes have no dynamical properties; e.g., they 
do not induce scattering! The remaining dynamics for modes below Aiow embody the effects 
of the high-momentum modes indirectly in the general spirit of the Renormalization Group 
method. 

It is evident that s has units of m~^; so we set sq = Aq^ and siow = Aj'^. In what 
follows, we take the MT-III potential for Vq and Sq and siow were chosen to coincide with 
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8 6 4 2 




FIG. 3: (Color online). The MT-III potential (x m/H^) in fm before (Vb(p, p'), upper contour) and 
after {Viow{p,p'), lower contour) SRG evolution from Aq = 16 fm^""^ to Aiow = 2 fm""'^. 



Ao = 16 fm~^ and Aiow = 2 fm~^, respectively. We obtain Viow by numerically integrating 
Eq. 21 from sq to siow using a 4tli order Runge-Kutta method or by using the Cayley 
approximation for the infinitesimal similarity transformation. Both methods give equivalent 
results. 3D plots of Vq and Viow are displayed in Fig. |3] where Viow has substantial off- 
diagonal contributions only below Aiow For larger momentum the off-diagonal components 
are suppressed and there is significant strength only along the diagonal which renders the 
potential sterile to scattering at larger momenta. These results confirm the general behavior 
reported by Bogner, Furnstahl, and Perry [17j. At low momentum Viow is substantially 
deeper than Vq. 

Next we perform a separable expansion of Viow The differences between the separable 
expansion of Vq and Vow are readily apparent in the Harms form factors of Fig. [2] where 
we show the first two form factors arising from the solution of Eq. |5] with V — > Vow- 
Fig. [2] also shows the same quantities for — > Vq, and reveals that the SRG suppresses 
the form factors of Vow at momenta larger than A^o^,. Moreover, there are no longer any 
form factors with large negative eigenvalues. The form factor associated with the bound 
state always has eigenvalue Ai = 1, of course, while the second has eigenvalue for Vow, ^2 = 
0.119368 . . . , is considerably different from the large negative value of A2 = —1.89166 . . . 
obtained with Vq. Also, with Vow, all remaining eigenvalues have magnitudes less than A2, 
e.g. A3 = 0.024633 . . . , and so on. (The SRG form factor associated with the negative 
eigenvalue, A = —0.11781 . . . , vanishes for q < Aq and is ignored.) This forms the basis for a 
systematic organization of terms in what can be considered a rapidly converging perturbation 
expansion [HI [19] . 

The calculations of the 2-body pcot 6 resulting from the rank-1 (UPA), rank-2, and rank- 
4 expansions of Vow appear in Fig. [T] along with the same results for Vq. We see that the 
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separable expansion of Viow for momenta befow Aiow = 2 fm~^ is essentially converged by 
rank-2 which is a modest improvement over the equivalent expansion for Vq while, by the 
nature of the similarity transformation, leaving the deuteron binding energy unchanged. 
These findings encourage us to employ the equivalent of the rank-2 T-matrix for View in the 
3-body, nd, elastic scattering calculations considered in the next section. As that discussion 
will indicate, for 3-body applications, we must be able to analytically continue the 2-body 
T-matrix away from the real axis without introducing spurious poles. We accomplish this 
by fitting the relevant form factors shown in Fig. [2] with sums of even polynomials of p times 
Gaussians. The fits we have obtained are excellent and would be indistinguishable from the 
numerical results in Fig. [2j 



IV. UNITARY POLE APPROXIMATION WITH A DIPOLE FORM FACTOR 

Before developing the higher rank separable expressions and to provide a model with 
which to compare, we briefiy review the example of a rank-1 separable 2-body s-wave treat- 
ment using a phenomenological dipole form factor closely following our earlier work [? ]. 
We limit our development to s-waves, but the forms easily generalize to the higher partial 



waves. Recalling Eq. 13 



V^-\^b){^b\ (22) 
where = 1 has been used. Now using a Yamaguchi or dipole form factor [33j, we have 

(p|V^B)=ArVVl>)=(^^^ (23) 

where g{p) is the form factor which satisfies g{0) = 1, p = \p\, and A/" is found by normalizing 
according to Eq. |6j 

= -Wb|G.(-B)Wb> = / tI^T^^ = 1^7^ (24) 

J [ZTiy B + q'^/m Svr (p + 7)"' 

where 7 = y/niB has been used. As discussed in Ref . [? ] , the dipole form factor is intimately 
related to the Hulthen wavefunction [3l]. Similarly, we find for the deuteron propagator. 



A-i(E+) = l-^^^±^^L_. (25) 
[3 + V-mE+ 



Finally, using the dipole form factor in Eq. [7] gives the wave function normalization factor. 
Using Eqs. 14, 23, and 25, we construct the fully off-shell elastic 2-body amplitude, /2, 

via 

Up,p';E) = T{p,p';E) ^ ~(^-^^gip') g{p")AiE)^ 



2 (^±11! I (l il+j)^) ' .27) 



/3 + V-mE+^ 
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and the on-shell scattering amplitude is given by 



{-PM2P + 7) + (3/3^ + 2/^7 + 7^)p^ + p^) 



-ip 



-1 



(28) 



Note that the 2-body amphtude has a pole at ii^ = — 7^/m as it must. We will see later 
that this pole structure plays a critical role in 3-body calculations by tracking the breakup 
threshold. In standard fashion we introduce the phase shift, 5, via 



f{p) 



p cot 5 — ip 



and see that Eq. [28] contains the correct unitarity term and that we can identify 



pcot 5 



+ , 3/3^ + 2/37 + f 2^ 
H ^^.^ — — p + 



2(/5 + 7)2 2/?(/3 + 7)2 
Comparison with the Effective Range Expansion, 

1 1 



;P'. 



pcot S 



allows us to identify 

2(/3 + 7)' 



/57(2/5 + 7)^ 



3/3^ + 2/^7 + 7^ 
/5(/3 + 7)2 



and Vi 



(29) 



(30) 



(31) 



(32) 



Note that the ERE terminates at 0{p'^) for the present choice of g{p). The exact 
deuteron binding energy for the MT-III potential is B = 2.23995 . . . MeV which implies 
7 = 0.232405 . . . fm~^. The numerical solution of the exact Lippmann-Schwinger Equation 
for this potential yields phase shifts which in turns can be used to compute p cot 5, a fit to 
which gives the values: 



,fit 



5.505 fm, r^^ = 1.891 fm and vf = 0.110 fml (33) 
Fixing 7 to the value given above and determining f3 = 1.281 fm~^ by fitting to a^* then 

0.170 fm^ 



Eq. [32] predicts 

ro = 1.897 fm and = 0.170 fm^ (34) 

which is in exceptional agreement with the fitted values considering the simplicity of the 
approximation. Similar levels of agreement are found using Gaussian form factors, g{p) = 
exp(— p^//5^), but the finite range expansion does not terminate. Note that the value of /3 
is roughly two pion masses which we would expect qualitatively for this "cutoff' parameter; 
i.e., (3 is "natural". Finally, we observe that none of the quantities in Eq. |32]is ill-behaved 
in the short-range limit {(3 00), namely, a — > I/7, tq ~ I3~^ 
expected. 



0, and f 4 ~ /3 



as 



V. FADDEEV EQUATIONS FOR "5^3 nd 

2 

Having established our method for treating 
problem of nd ^Sz (or "quartet") scattering. 



SCATTERING 

the 2-body interaction, we now address the 
In an earlier work [? ], we made use of 



12 



FIG. 4: (Color online). Diagrammatic representation of the Born ("pinball") term of 3-body 
particle-dimer elastic scattering Faddeev equation. The " blobs" are the form factors which depend 
on the magnitudes of the relative 2-body momenta which are |p/2 + q| and |p + q/2| for the lower 
and upper form factors, respectively. 



the Unitary Pole Approximation, rank-1 separable approximation, to study convergence 
properties and limitations in treating finite range effects using effective range expansions 
(ERE's). In that work we found that ERE's can account for finite range effects adequately to 
Next-to-Leading (NLO) order, but fundamentally different terms, similar to 3-body contact 
terms, arise in the finite range expansion of the UPA at Next-to-Next-Leading order (N^LO) 
that do not appear in the ERE formulations. Furthermore, in treating 3-body systems, there 
are catastrophic failures of unitarity in the form of negative imaginary phase shifts in low 
order ERE treatments. Including vertex corrections from explicit pions treated in chiral 
perturbation theory have been shown to cure this disease [35], but at a consider cost in 
effort and transparency. 

We now construct the 3-body Faddeev equation in the Unitary Pole Approximation. 
Consider the scattering of a low energy neutron from a deuteron in the quartet, ^Sz , channel. 
Treating the nucleons as identical particles in an isospin formalism, the S-wave (L = 0) 3- 
body t-matrix is given by [22] 

E) = 2Z^^^ip,p'; E) - j q'dq 2Z^^\p, q; E) ^{E - ^) X^^^(g,p'; E), 

(35) 



where ^{E) is given by Eq. 25, identified with the deuteron (or dimer) propagator, and Z 



is the particle-exchange ( "pinball" ) amplitude depicted diagrammatically in Fig. |4j, 

Z-^^M, E) ^ -A, - /' ,, VB(|P + P72|)fa(|p' + P/2|) ^ 

2 J-1 P + P ~ TnE+ + pp X 

where Pe{x) is a Legendre polynomial and Ajs is the isospin-spin structure factor which is 
unity for bosonic scattering. For the nd quartet case, 

A.^.Q,l,|i;M)..(n,n;O,o).0x(-i) (3, 

where the U^s are the unitary Racah recoupling coefficients. Examining the homogeneous 
terms, we note that by associating the iPb{p) factors from Z{p, q; E) with the dimer propa- 
gator, A(i^^ — ^), the (off-shell) 2-body T-matrix, Eq. 14, is reconstituted. 
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FIG. 5: (Color online). The real part of L = quartet phase shifts using three implementations 
of the rank-1 (UPA) model. The solid curve (black) is the rank-1 calculation for Viow and the dashed 
curve (blue) is the rank-1 calculation for Vq. The dotted curve (red) is the rank-1 calculation using 
the phenomenological dipole form factors discussed in Sec. |IV| Also shown is the effect of reducing 
the cut-off reduced from A = 16 fm~^ to A = 2 fm~^. This change has no significant effect on 
the rank-1 calculation for Vq and yields the dash-dot-dot curve (blue) is the rank-1 calculation for 
Viow The triangles are the full Faddeev calculations of Huber, et al. |12J and Kievsky, et al. [T3] 
and references therein. 

Including the wavefunction renormalization factor, Eq. |26} defines the off-shell 3-body 
T-matrix, 

T^^ip^p'; E) = ZBX^^^ip,p'; E), (38) 
The Faddeev equation for T3 becomes, 

T^''^{p,p'-E) = Ti--{p,p'-E)+ / dqq'K^''^{p,q-E)T^^^{q,p'-E), (39) 

^0 
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FIG. 6: (Color online). The imaginary part of L = quartet phase shifts using three imple- 
mentations of the rank-1 (UPA) model. The solid curve (black) is the rank-1 calculation for Vio^ 
and the dashed curve (blue) is the rank-1 calculation for Vq. The dotted curve (red) is the rank-1 
calculation using the phenomenological dipole form factors discussed in Sec. |IV[ Also shown is the 
effect of reducing the cut-off reduced from A = 16 fm~^ to A = 2 fm~^. This change gives the 
dash-dot curve (black) for the rank-1 calculation for Vq and the dash-dot-dot curve (blue) for the 
rank-1 calculation for Vq. The triangles are the full Faddeev calculations of Huber, et al. |12j . The 
breakup threshold is at p = 53 MeV/c. 



where 

r3^°™(p,p';E) = 2ZBZip,p';E), and (40) 



dx i)B \ \ + -r + pqx ]i'B\ A/ "T + <? + PQ^ ^UPA E 



J_i ^ \\j 4 / \V4 / V 4m 

(41) 



Using the dipole form factor, Eq. 23 , we can perform the needed angle integration in Eq. 
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to find the UPA result, 



Ti'^^{p,q-E) 



IS 



3mpq{p^ - g2)(3g2/4 -(3^- mE){3p^/A ~ (3'^ - niE) 
3 



-{p — g ) In 



p^ + q^ — mE + pq 



- {3p'^/A- (3^ -mE)\n 
+ {3qy4- (3^ -mE)\n 



p2 _|_ g2 _ _ pq 

+ + 4/5^ + Apq 



4p2 + g2 _^ 4^2 _ 

p2 + 4g2 + 4/^2 _^ 4pg 



p2 _^ 4^2 _^ 4^2 _ 



(42) 



Likewise for tlie 3-body kernel, we find 

4/?(7 + /3)2A, 



Snmpq 



{l3+y/3q'2/4-mE)^ 

3 



(3p2/4 - /52 - m£;)(3gV4 - /^^ - mE) 



In 



p2 _^ g,2 



mE + 



p2 _|_ g2 _ _ pg 



1 

-g2)(3p2/4_^2 



In 



In 



4p^ + g^ + 4/?^ + 4pg 



4p2 + g2 + 4^2 _ ^pq 
p^ + 4g2 + 4/3^ + 4pg 



p2 _^ 4g2 _^ 4^2 _ 4pg 



The on-shell 3-body S-wave scattering amplitude in UPA is given by, 



UPA 



m 



on 



3p2 
4m 



11' 
m ' 



(43) 



(44) 



In Ref. [? ], this result is used to study the convergence of finite range expansions and 
compare with similar expansions using effective field theories. Here we extend this to higher 
order in the separable expansion. 



VI. BEYOND THE UPA 



To move beyond the UPA we examine the 3-body scattering equations for separable 
potentials of arbitrary, but finite, rank as introduced in the 2-body context in Section |TT} 
We closely follow the treatment of Watson and Nuttall with minor changes in notation. We 
restrict our attention to the scattering of one particle from the bound state ("dimer") of 
the other two. We work in the center-of-mass system and assume that the particles have 
identical mass, m. To begin, we recognize that the 2-body form factors now act in the 
3-body Hilbert space; so a more elaborated notation is needed that allows us to identify 
the various 2-body partitions. We adopt the "odd-man-out" notation where the subscript i 
refers to the {j7c}-pair where neither j nor k equals z, e.g. 3 refers to the {12}-pair. Thus, 
the 3-body state with particles 1 and 2 in state \ipn) and particle 3 in state \p^) will be 
represented by 

|7/;„(1,2))® Ips) ^ |n3,p3). (45) 
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FIG. 7: (Color online). The real part of L = quartet phase shifts using three implementations 
of the rank-2 Harms separable expansion. The solid curve (black) is the rank-2 calculation for Vq 
and the (indistinguishable) dashed curve (blue) is the rank-2 calculation for Viow The dotted curve 
(red) is the rank-2 calculation using phenomenological (gaussian times polynomial) form factors 
fit to the 2-body phase shifts. For this case the effect of reducing the cutoff from A = 16 fm~^ to 
A = 2 fm~^ has no significant effect. The triangles are the full Faddeev calculations of Huber, et 
al. [12] and the crosses are the calculations of Kievsky, et al. [13j (and references therein). 

Let E be the 3-body energy and G''q\e) be the free 3-body propagator, then 



where Gq (z) is the free 2-body propagator and /ij is the reduced mass of particle i with the 
7-dimer. Under the assumption of equal masses, fii = 2m/3 independent of i. If B is the 
dimer binding energy, then the particles are on the energy shell when E = Sp'^/Am — B and 



{n^,p,'\Gi^\E)\nj,Pi) = {2nr6^'\p,' - p^){n^\G^^\E - /2^^,)\nJ) , 



(46) 



{n',\G^^\E-3py4m)\n,) ^ {n',\G^^\-B)\n,) = 



(47) 
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FIG. 8: (Color online). The imaginary part of L = quartet phase shifts using three implemen- 
tations of the rank-2 Harms separable expansion. The solid curve (black) is the rank-2 calculation 
using for Vq and the dashed curve (blue) is the rank-2 calculation for Viow The dotted curve (red) 
is the rank-2 calculation using phenomenological (gaussian times polynomial) form factors fit to 
the 2-body physics. Also shown is the effect of reducing the cutoff from A = 16 fm~^ to A = 2 
fm^^. This change has no effect on the Viow calculation, but gives the dash-dot (black) curve for 
Vo which violates unitarity below the breakup threshold at p = 53 MeV/c. The triangles are the 
full Faddeev calculations of Huber, et al. [12] and the crosses are the calculations of Kievsky, et 
al. [13] (and references therein). 
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FIG. 9: (Color online). The real part of L = 1 rwi quartet phase shifts using three implementations 
of the rank-2 Harms separable expansion. The solid curve (black) is the rank-2 calculation using 
for Vb and the dashed curve (blue) is the rank-2 calculation for Viow The dotted curve (red) is 
the rank-2 calculation using phenomenological (gaussian times polynomial) form factors fit to the 
2-body phase shifts. The triangles are the L = 1 (^-Pi/2, '^P3/2^ ^-^5/2) full Faddeev calculations of 
Huber, et al. [12] and the crosses are the calculations of Kievsky, et al. [13] (and references therein). 



as in Eq. [6]of Section [TTj The Faddeev Equation for the 3-body wave function for, e.g., 
is (see e.g. Watson and Nuttall, Eq. 4.40) 



/|<!p3;.)\ 




Ti{E) TiiEy 

mE) mE) 

T-s{E) T-,{E) 



',(1) 



(48) 
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200 



FIG. 10: (Color online). Same as for Fig. [ojbut for the imaginary part of the phase shifts. The 
breakup threshold is at p = 53 MeV/c. 



where in {ip^^ „ .^), the z-pair interact last so the i-th particle is free asymptotically with 



momentum pa and where Tj{E) are 2-body T-matrices for the z-pair acting in the 3-body 
space which for a rank-N separable potential are given by, 



T,{E) = - E Kp'!) A„,„,(E - 3pf /4m) (m^p^'l 



(49) 



where Anmiz) appears in Eq. 11, Hence 

TV 



<!p3;i.) = -G^o\e) l^l'Pi) A,„,(i5;-3pf /4m) 

li,mi 



X 



("^l,P;i<!p3;i.) + K>Pll<!p3;i.) 



(50) 
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(2) X 



= -G^^\e) E / ^ 1/2, P'^') A.,„,(^-3pf/4m) 



X 



(^2,P'2VSp3;s) + (^2,P'2V: 



and 



X 



(^3,p'3l<!p3;i.) + (^3,p'3'l<!p3;^) 



Now define the 3-body amplitudes, 



to find explicitly, 



-'^n'=n3(P3'P3;-E') 



-E 



(51) 



(52) 



(53) 



(2,m2 



^ (n'3,p^|G(^)(£;)|Zi,p;') A,,^,(£;-3pf/4m) ^^.^^(p;', P3; £;) 

^ {n',,p',\G^^\E)\k,p'^ A,,^,(£;-3pf/4m) ^^.^^(p',', P3; £;), 

(54) 



and 

-'^<n3(P2>P3;-^) 



X^^^.^{p[,Ps;E) = (n'i,p;|Gj')(£;)|n3,P3) 



^^'i'Pil^o'Hi^)K2,p'2') A,,„,(E-3pf/4m) ^^^.^(p;', P3; S) 



l2,m2 



E 



^V3' / , 



{n'-„p[\G^^\E)\k,p'^) Aum-SE-2,pf/Am) ^^(p'g', pa; £:), 



//2 



(27r)= 



(55) 



= (n'3,p^|4=^)(£;)|n3,P3) 



•E 



{n'-^,p'^\Gf{E)\h,p'[) Ai-^m,{E-3p'('/4m) X„(p;',p3;^) 



//2 



(27r)^ 



{n'2,p'2\G^o \E)\k,p'^) Ai,m,{E-?,pf/Am) X^3n.(P3, P3; £^)- 



(56) 
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So far, the three particles have been treated as distinguishable. Now we assume they are 
identical spinless bosons in which case we may, without loss of generality, let 



X^. (p:, p, ; E) ^ 5.,X^„(p', p; E) 



(57) 



and 



X„,„.(p:,p,;^), X„.,„.(p:,p,;E) ^ (1 - 5,,-)^Jn(p', P; ^) (58) 
where and X-^ are "direct" and "exchange" amplitudes, respectively. We then have 

X^'niP'.P-.E) = -2f2 1 1^ {n',p'\G^^\E)\l,p")x - ?>p"' / Am) X^„(p",p;E), 

(59) 



and 



X^„(p',p;E) = {n\Y>'\Gt\E)\n,Y>)x 

N 

-E 



Lm 



{n', p'\Gf\E)\l, p")x AUE - 3p"'/4m) 



(27r) 



where the "X" subscript in (n', p' |G(=^)(i^)l^,P)x means that the bra and ket must corre- 



X 



^mnip"^ P; E) + X^„(p", p; E) 



(60) 



spond to different arrangements of the three bosons as required by Eq. 54 56 Then define 
X = X^ + 2X^, whence the two previous equations can be combined to yield 



X„,„(p', p; E) = 2{n', p'\Gi^\E) \n, p)x 



(61) 



-2E 



Lm 



d^p" 

(27r)3 



{n',p'\G'i\E)\l,p")x AUE-3p"'/Am) X„„(p",p,E). 



//2 



We now write 



n',p'\G^^\E)\n,p)x 



(^n'|p^ + f)(p+^|^n) 

e-Up'+p" + p'-p) 



(62) 



and then define partial wave components, 

"+1 



Z^,^{p\p;E) = ^ dxPLix) {n',p'\G^^\E)\n,p) 



X 



E — —{p + + p'px) 



(63) 



where Pl{x) is a Legendre polynomial. There are then, for a rank-A^ separable potential, 
A^^ coupled Faddeev equations for bosonic scattering of the form. 



X„\(p',p;E) = 2Z^,„(p',p;E) 
2 ^ 



27r2 



poo 

J2 / P"'dp" Z}:,,{p\p"; E) AUE - 3p"'/Am) X^p'^p; E). 

Lm 



(64) 
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The corresponding equations for nd quartet scattering are obtained by the simple expedient 
of multiplying the Z's by the isospin-spin structure factor A/5 = —1/2 (see Eqns. 36 and 37). 

We now consider the special case of a rank-2 separable potential, and we observe that 
the 2^ equations decouple into two pairs of coupled equations. We will focus on elastic 
scattering, so, assuming that n = 1 corresponds to the deuteron ground state, we focus on 
the equations for Xu and X21: 



X'^,{p\p-E)= 2AjsZ^,{p',p;E) 



2A/5 



27r2 



1/2 7 // 

p dp 



Z^^{p\p";E) An(^-3/74m) 



+Z^^{p\p"- E) A,,iE - 3/74m) 



Xt,ip",p;E) 



2A 



IS 



27r2 



p"^dp" 



Z^^{p\p"-E) Ai2(^ - 3p"74m) 



+Z'^^{p\p"- E) A22{E - 3/V4m) 



X^,{p",p;El 



(65) 



and 



X^,{p',p;E)= 2AjsZ^,ip',p;E) 



2A 



IS 



27r2 



p"^dp" 



Z^,{p',p";E) AniE-Sf ^Am) 



+Z^,ip',p";E) A2i(i?-3/74m) 



Xi^,ip",p;E) 



2A 



IS 



27r2 



1/2 7 // 

p dp 



Z^,{p'y-E) A,2{E-3p"'/Am) 



+ZUp'y: E) A22{E - y/Am) 



Il2 



X^,{p\p;E). 



(66) 



The nd quartet elastic scattering amplitude for the L-th partial wave is 



3p2 



Am 671 Am 



3p2 



(67) 



where Zi = Zb is the bound state normalization given by Eq. [7j 



VII. RESULTS 

We begin with calculations of the ^Ss phase shifts using three implementations of the 
rank-1, unitary pole approximation (UPA). Shown in Figs, [s] and [6] are calculations of the 
real and imaginary ^Ss phase shifts for our three implementations of the UPA compared 
to the full Faddeev calculations of Huber, et al. |12] and Kievsky, et al. [13j and references 
therein. The three implementations are (1) the phenomenological dipole form factor treated 
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in Section IV, (2) the rank-1 Harms form factor for the bare Malfleit-Tjon triplet potential 
(Vq), and (3) the rank-1 form factor for Viow which is obtained through a similarity renormal- 



ization group evolution of Vq from 16 fm ^ to 2 fm ^ as described in Section III All three 



implementations accurately reproduce the real part of the s-wave quartet phase shift of the 
full Faddeev calculation, and each method gives a respectable reproduction of the imaginary 
S-wave phase shifts when the cutoff is large, Aq = 16 fm~^ . Also shown is the effect of 
reducing the cutoff momentum in the 3-body calculation from Aq = 16 fm~^ to Aq = 2 
fm~^. The real part of the phase shift is relatively insensitive to this change. Reducing the 
cut-off momentum causes the imaginary part of the phase shift to be reduced somewhat at 
the larger momenta, and, in the case of Vq introduces a spurious negative imaginary phase 
shift at low momenta. No such spurious effect is seen using Viow however. 

Shown in Figs. Mand 8 are calculations of the real and imaginary ^Ss phase shifts for three 
implementations of rank^ separable expansion compared to the full Faddeev calculations of 
Huber, et al. [I2] and Kievsky, et al. [13] and references therein. The three implementations 
are (1) a phenomenological form factor (polynomial times gaussian) fit to the NN phase 
shifts, (2) the rank-2 Harms expansion using the first two form factors for the bare Malfleit- 
Tjon triplet potential (Vq) and (3) the rank-2 Harms expansion using the first two form 
factors for Viow which is obtained through a similarity renormalization group evolution of Vq 
from 16 fm~^ to 2 fm~^ as described in Section HI The form factors are shown in Fig. |2} 
Each method gives a respectable reproduction of the full Faddeev calculations for both the 
real and imaginary S-wave phase shifts. Also shown is the effect of reducing the cut-off 
momentum in the 3-body calculation from Aq = 16 fm~^ to Aq = 2 fm~^. Similar to the 
rank-1 case, we see that reducing the cut-off momentum does not significantly affect the real 
part of the phase shift and causes the imaginary part of the phase shift to be reduced at the 
larger momenta but less so than in the rank-1 case, and, as in the rank-1 case, introduces 
a small spurious negative imaginary phase shift at low momenta for the Vq implementation. 
We conclude that the high momentum contributions are necessary to retain unitarity below 
threshold for the Vq case, but these contributions are included in Vow such that unitarity is 
respected even when the cut-off is reduced. 

We have also calculated higher partial waves and find similar levels of agreement and 
systematics. Shown in Figs, [o] and 10 are calculations of the real and imaginary ^Pj phase 
shifts for three implementations of rank-2 separable expansion compared to the full Faddeev 
calculations of Huber, et al. [12] and Kievsky, et al. [13] and references therein. Since our 
simple model does not have spin-orbit coupling, our phase shifts are independent of the j- 
value. Nevertheless, the calculations of the real part of the phase shift lies within the range of 
results from the full Faddeev calculations including spin-orbit forces while the results for the 
small imaginary part of the phase shift fall just below those of the full Faddeev calculations 
at the higher momenta. 

In Table [TT] we present calculations of the total, elastic, break-up, and back-angle cross 
sections for quartet nd scattering at Eiab = 10 MeV summing partial wave contributions 
to Lmax = 8 compared with the results of Koike and Yaniguchi [36]. The back-angle nd 
cross section is dominated by the quartet term. As already observed long ago, for these 
quartet-dominated scattering measurements, the UFA results are quite respectable, espe- 
cially considering their simplicity. The quality and similarity of the three rank-1 calculations 
is somewhat surprising considering how different are the various rank-1 form factors, shown 
in Fig. [2| Apparently, including the correct 2-body phenomenology through the leading form 
factor accounts for most of the physics involved in quartet scattering. Still, it is satisfying 



24 



Model 


atot (mb) 


aei (mb) 


asR (mb) 


^{9 = 180°) (mb) 


UPA (Form Factor) 


752.3 (751.6) 


712.9 (712.2) 


39.4 (39.4) 


177.6 (177.4) 


UFA {Vo) 


776.4 (777.7) 


736.9 (739.6) 


39.5 (38.1) 


179.8 (180.5) 


UFA (Viou.) 


789.6 (789.5) 


749.6 (749.6) 


40.0 (40.0) 


182.0 (182.0) 


Rank-2 (Gaussian x polynomial) 


770.9 (770.9) 


730.9 (731.3) 


40.0 (39.6) 


179.2 (179.4) 


Rank-2 Harms (Vq) 


779.0 (780.3) 


739.3 (742.1) 


39.6 (38.2) 


180.9 (181.6) 


Rank-2 Harms (Viow) 


794.4 (794.3) 


754.3 (754.2) 


40.1 (40.1) 


184.2 (184.2) 


Koike &Yaniguchi |36] 


812.3 


770.2 


42.1 


185.0 



TABLE 11: Calculations of center of mass total, elastic, break-up, and back-angle cross sections 
for quartet nd scattering at Eiab = 10 MeV, summing partial waves to L^ax = 8 using a cutoff of 
Ao = 16 fm~^. (The results from reducing the cutoff to Aq = 2 fm~^ are given in parentheses.) 



that including the rank-2 term systematically improves the results for all cases considered 
and the rank-2 Harms expansion of the SRG-evolved potential, Viow, yields the best results 
compared with full Faddeev calculations. We also note that the cross sections calculated 
using the rank-2 Harms expansion for Viow are largely insensitive to the cutoff, giving essen- 
tially identical results down to Aq = 2 fm~^. 

VIII. CONCLUSIONS 

In this work we reviewed Harms' method of separable expansions and, for the first time, 
applied it to a potential evolved to low momentum using the similarity renormalization 
group. This work differs significantly from that of Kamada, et al. (who use Legendre 
polynomials for the separable form factors of an RG-evolved potential) in that our Harms 
form factors, derived from the 2-body potential, incorporate the 2-body physics, and in 
particular the bound state physics, in the leading term. The goal of this work was to develop 
a computationally simple yet accurate approach to 3-body calculations. To demonstrate the 
utility and feasibility of this approach while avoiding the complications of tensor coupling and 
Efimov physics, we applied the Harm's separable expansion to 3-body nd quartet scattering 
using both the bare Malfiiet-Tjon-III potential, Vq, and its SRG-evolved potential, Viow We 
find that excellent agreement with full Faddeev calculations can be achieved at just rank-2 
and, not surprisingly, that the results using Vow are stable with respect to the momentum 
cutoff avoiding spurious non-physical negative imaginary phase shifts that arise for Vq when 
the cutoff is lowered. 

The nd quartet case is an important proof-of-principle example but has has long been 
known to be well described by the unitary pole approximation, so the quality of our rank-1 
results and improvements from higher rank terms may not seem surprising. The real test of 
our approach will come in future applications to the nd doublet case which is theoretically 
more challenging due to complications of tensor coupling and Efimov physics. 
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